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Abstract 

The bb production cross section in 920 GeV proton-nucleus fixed target collisions is measured by observing double 
muonic decays of ^-flavoured hadrons in the kinematic region —0.3 < Xf(^i) < 0.15. A total number of 76 ± 12 bb 
events is obtained with a likelihood fit of the signal and background simulated events to the data. The resulting cross 
section is a b i = 16.2 ± 2.5 sta t ± 2.8 sys nb/nucleon, or, when combined with a previous HERA-B measurement 
of similar precision, a b i = 15.4 ± 1.7 sta t ± 1.2™ corr - ± 1.9 c s ° r s r - nb/nuclcon, which is consistent with recent NLO 
calculations. 
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1 Introduction 

The measurement of bottom production in fixed target col- 
lisions offers the possibility to test perturbative QCD in the 
near threshold energy regime, where the effect of higher 
order processes, such as soft gluon emission, has been cal- 
culated CD El- At first order, the production mechanism at 
the HERA-B energy (y/s = 41.6 GeV) is dominated by 
gluon-gluon fusion (gg — » bb) 01. 

Three experimental results are published. Two are in- 
consistent, even though they were obtained by similar ex- 
periments searching for J/^ and semi-leptonic [5] de- 
cays of 6-flavoured hadrons. HERA-B recently published 
the most accurate result based on a measurement of J ftp 
decays of the ^-flavoured hadrons @. 

In this paper, a measurement of the bb production cross 
section performed with a b tagging technique independent 
from our previous measurement is presented. 

After production, bb pairs hadronise and mostly decay 
into c-flavoured hadrons. Since ^-flavoured and c-flavoured 
hadrons have a large probability to decay with the emission 
of a muon ("semi-muonic decay") Q, the bb production 
cross section is measured by searching for bb — > \i\x + X 
decay events, in which at least two of the four heavy quarks 
typically produced in a bb event (b, b, c, c) undergo semi- 
muonic decays ("double muonic b decays"). 

The bb event selection is based on a pair of oppositely 
charged muons not coming from the primary interaction 
vertex, having a large momentum transverse to the beam. 

2 Detector and data sample 

HERA-B (HUH is a large acceptance forward spectrometer 
installed at the 920 GeV proton storage ring of DESY. The 
Feynman-x (xp) of accepted muons from semi-muonic b 
decays is between —0.3 and 0.15. The detector is used to 
reconstruct charged particle tracks produced in the interac- 
tions of the proton beam halo with wires of different mate- 
rials ( 12 C, 48 Ti and 184 W), in several configurations ifTOl . 
Particles are tracked with a silicon microstrip detector iPTTTl 
whose first station (of 8) is a few centimeters from the tar- 
get system and which extends approximately 2 m further 
downstream. A primary vertex resolution of 500 /xm along 
the beam and 50 fim in the perpendicular plane is achieved. 
Up to 13 m downstream of the target, honeycomb chambers 



in the outer region lfl"2l[T3l . and microstrip gaseous cham- 
bers in the inner region lfi"4l . allow to track particles and to 
measure their momenta from the bending in a 2. 13 T-m ver- 
tical magnetic field. A Cherenkov detector lfl5l is used for 
n/K/p separation. An electromagnetic shashlik calorime- 
ter lPT6l serves for e and 7 identification. At the rear of 
the detector, muons with momenta larger than 5 GeV/c are 
tracked with triple stereo layers of gaseous tube chambers 
interleaved with hadron absorbers ifTTl . 

Double muonic b decays are searched for in 164 • 10 6 
events of p-C, p-Ti and p-W interactions (1.4 overlapped 
interactions per event, in average) collected with a mul- 
tilevel dilepton trigger HH in the 2002-2003 data taking 
period. The trigger is designed to select dilepton decays of 
J/tp mesons produced in the proton-nucleus collisions. By 
applying dimuon selection criteria similar to those of ref- 
erence [6], about 146,000 prompt J/ip mesons are recon- 
structed. With this event selection, double muonic b decays 
are a tiny fraction of the surviving muon pairs (Figured]). 
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Figure 1. Invariant mass of opposite-sign muons after dimuon 
selection (white histogram). Two peaks emerge above the back- 
ground (J/tp and ip'). The fit of the 3/tp signal at 3.097 GeV/c 2 
(grey line) includes events in the radiative tail [ 19 1. The grey his- 
togram is a Monte Carlo simulation of double muonic b decays 
(see text for details). 

The systematic uncertainties due to the detector per- 
formance and acceptance is reduced by normalising the bb 
production cross section to the prompt J/ip cross section. 
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3 Monte Carlo simulation 

The acceptance and the reconstruction efficiencies are de- 
termined through Monte Carlo simulations of the physics 
processes occurring in proton-nucleus interactions. For bot- 
tom and charm production, PYTHIA 5.7 |[20l is used. For 
bottom production, the predictions of NRQCD models are 
used, while 3 ftp production is tuned to match the xf and 
Pt distributions measured by other experiments. The frag- 
mentation process is simulated by JETSET 7.4 [20]. The 
energy left from the hard scattering is used by FRITIOF 
11211 to simulate the underlying inelastic event. The re- 
sponse of the detector is simulated by GEANT 3.21 [22]. 
See, for more details, reference (6l. 

4 Measurement method 

The bb production cross section can be expressed as 



N± Affj/^BRj/^ 



bb 



X 



bb 



where is the total number of double muonic b decays 
in our sample, Actj/^, = 417 ± 37 nb/nucleon ll23l 

is the prompt J/tp production cross section in the detector 
acceptance (—0.35 < xf(J/"0) < 0.15) and BRj/^, is the 
branching ratio for J/ip — > decays (5.88 ± 0.10%) 

iPTl . The term is defined as 



X 



bb 



£,BR JJi (l-0 i )e 



bb,j,i 



% 



where the index i runs over 14 different target configura- 
tions. The index j refers to the different ways to generate 
dimuons via semileptonic b or b — > c decays (Section [5]) 
having different branching ratios (BR#j) and reconstruc- 
tion efficiency (£^~ { , the superscript indicates the dimuon 
charge). The charge factors 9j gives the fraction of result- 
ing same-sign muons (Table [D, Nj M{ is the number of 
J/ip mesons reconstructed with efficiency eju^, A is the 
atomic weight of the target and a = 0.96 ± 0.01 [24J is the 
J/ip nuclear suppression in the central xf region. Nuclear 
effects are expected to be negligible for open bottom and 
open charm production. 



The number NX is obtained with two methods lead- 

bb _ 

ing to two determinations of the bb production cross section 
(Section El). 



5 Signal decay modes 

The signal sample consists of events with two muons com- 
ing from semi-muonic decays of heavy quarks. Four heavy 
quarks are typically produced in a bb event (b, b, c, c). 
Depending on the type of hadrons decaying in the semi- 
muonic mode (bb, bc+bc, bc+bc or cc), four classes (Bj) of 
signal events are defined ll25l (Table [j}- 



Class 


Decaying hadrons 


Branching ratio 


6 


Bl 


bb 


0.0084 


0.187 


B2 


bc+bc 


0.0160 





B3 


bc+bc 


0.0156 


0.742 


BA 


cc 


0.0081 


0.315 



Table 1. Branching ratios and charge factors 9 (see text for de- 
tails) of four classes of double muonic b decays. The branching 
ratios are correlated and are affected by a relative uncertainty of 
10%. The relative uncertainty on is less than 1%. 

For each class, the branching ratio for dimuon decays is 
obtained from the branching ratios for semi-muonic b and 
c decays reported in the PDG [7 [. The fraction of decays 
into same-sign muons (/i ± /x ± ) is given by a charge factor 6, 
1-9 being the fraction of decays into opposite-sign muons 
(^ + /i~). The value of 9 is determined through Monte Carlo 
simulations and includes the effect of B and B$ mixing 
and all possible decay paths leading to muon pairs. 

A fifth class, which is not included above, must be also 
considered. Events in which pairs originate from 

double muonic decay of cc pairs from b — ► ccs + X de- 
cays represent 5.5% of the total signal events. These events 
are assigned to the class B2, which exhibits a similar fi- 
nal state. The systematic uncertainty corresponding to this 
choice is included in the branching ratio BR(6 — ► ccs + X). 

Events with more than two muons might fall into more 
than one class. In order to avoid multiple counting of signal 
events, the calculated branching ratios in Table Q] include 
the probability that 6-flavoured and c-flavoured hadrons do 
not decay into muons, and events are assigned to classes 
with a priority given by the order of classes in Table Q] 
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6 Background contributions 

The main background for 66 event selection is made of 
events having two muon-like tracks not coming from the 
primary interaction vertex. Two sources of this type are 
considered: double muonic decays of c-flavoured hadrons 
and random combinations of muons from decay of low 
mass mesons (mainly pions and kaons). The latter back- 
ground is referred to as "combinatorial background". 

The Monte Carlo simulation shows that a cut around 
the 3/i/j mass (between 2.95 GeV/c 2 and 3.25 GeV/c 2 ) re- 
moves background events from b — > J/ip + X decays and 
prompt J/tp decays. The 3/ip mass region is excluded to 
be also statistically independent from the measurement of 
bb production cross section in Q, where bb events are iden- 
tified with J/tp decays of the ^-flavoured hadrons. 

Muons from Drell-Yan events are at least 10 times less 
abundant than the combinatorial background in the invari- 
ant mass region of interest [26]. 

The number of background events from double muonic 
c decays (N^) is determined through a Monte Carlo sim- 
ulation normalised to the number of prompt J ftp mesons 
reconstructed in the data. It can be expressed as 



X 



+- 



Acrj/^BRj/^, 



-X c 



where a C c is the charm cross section [27] (49±5 /xb/nucleon 
(21), BR cg = BR(c fi + X) 2 = (0.082 ± 0.005) 2 0, 
and the term X c5 is defined as 



X r . 



3/ip,i 



A 



l-a 



'3/ip,i 



The term e c ^i is the bb selection efficiency for double 
muonic c decay events in the target configuration i. 

The combinatorial background in the channel 
(iV+~) is determined with /i ± /x ± data. 

Muons from double muonic b decays do not come from 
the same decay vertex. However, the forward boost in fixed 
target collisions is such that tracks coming from two long 
lived particle decays are almost as close as those origi- 
nating from a single particle decay. In double muonic b 
decays, the middle point of the segment of minimum dis- 
tance between the two muons (pmd) is preferentially lo- 
cated downstream of the target, while the region upstream 
of the target is dominated by combinatorial background, in 
both n + n~ and /U^/z^ channels. 



Assuming that the combinatorial background in the two 
final states has a similar shape, the number of events are 
normalised with respect to the upstream side and the com- 
binatorial background in the /U + /i~ channel is estimated as 
the difference between the number of // ± ^ ± pairs in data 
(iV~) and those expected from double muonic b decays 
in the same channel (Njj^) 



N, 



N, 



bb 



is obtained from Monte Carlo simulations, under 

bb 

the assumption that a b i = 15 nb/nucleon, with the formula 



N, 



CT; 



bb 



A^j/^BRj/v 



-X 



±± 

bb ' 



where the term X^ is defined as 



bb 



i_ a Dj BR gj g j £ S,i 



x bt = E N ^ A 



'3/ip,i 



The term er^ . indicates the reconstruction efficiency 
for double muonic b decays into // ± ^ ± . 



7 Data analysis 

Two methods have been used to estimate the number of 
events from bb decay. A first method, which is used to ob- 
tain our final results, is described in Section 1731 and a sec- 
ond method, which is used as a cross-check, is described 
in Section 1741 With one exception (discussed below), the 
choice of cuts used in both methods is the same. Initial 
muon selection cuts are given in Section 17.11 and the pro- 
cedure for optimizing the final cuts is described in Section 



The selection of 66 decay events begins by requiring 
that the events have at least two muons. 

7.1 Muon selection 

A first general muon selection is performed by requiring 
a high-quality reconstructed triggered track having a mo- 
mentum between 5 and 200 GeV/c, a minimum transverse 
momentum (pr) of 0.7 GeV/c and a minimum x 2 proba- 
bility of the track fit (P x ) of 0.003. The muon likelihood 
|[29l . as measured in the muon detector, must be greater 
than 0.05. 
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7.2 Cut optimisation 

Both muons of an event are required to have a minimum px 
of 1 GeV/c and a minimum impact parameter to the target 
(Ip) of 1.5<7 (where a is the Ip resolution). The Ip is de- 
fined as the perpendicular distance between the target wire 
and the point on the track extrapolated to the z -position of 
the target. The dimuon invariant mass (m) is required to be 
at least 2 GeV/c 2 . 

The optimal bb selection criteria are found by maximis- 
ing the signal significance S, which is defined as 



y/N±- + N+- + N+-' 

The number of signal events (N^~), which consists of 
oppositely charged muons, is obtained from Monte Carlo 
simulations, under the assumption that a b i = 15 nb/nucleon: 

N±~ = ^ X±~. 

The formulae used to estimate the number of back- 
ground events (N^~ 5 ~ and N+~) have been presented in 
Section [6] 

A large fraction of combinatorial background consists 
of muonic decays of kaons and pions. Since the angle be- 
tween the emitted muon and the decaying particle (kaon or 
pion) is small, such a background is suppressed by increas- 
ing the lower limit on P x . 

An upper limit on the kaon likelihood (L^) 11301 . as 
measured in the Cherenkov detector, suppresses muon can- 
didates from kaon decays. 

A lower limit on pt suppresses muons from low mass 
particle decays, since they are expected to have a smaller 
Pt than those from b decays. 

A lower limit on Ip discriminates muons originating 
in b decays from background muons. Since the Ip is corre- 
lated with the lifetime of the decaying particle, and the life- 
time of ^-flavoured hadrons is larger than that of c-flavoured 
hadrons, the impact parameter cut suppresses open charm 
background. 

In order to suppress background in the proximity of the 
target, a lower limit on the difference between the pmd and 
the target positions along the z-axis (Az) is applied. 

Unphysical events are suppressed by the requirement 
of the following upper limits: Ip < 50a, pr < 5 GeV/c, 



|Az| < 5 cm and m < 8 GeV/c 2 . The optimisation of the 
last three cuts (lower limits on pr, Ip and Az) is performed 
simultaneously. 

The optimisation procedure for S results in the bb se- 
lection criteria listed in Table [2 where the number of sur- 
viving and fi^^ pairs in the data, at each selection 
step, are also shown. After applying all cuts, the number of 
remaining dimuon candidates is 

AT+- = 167±13 stat . 



Cut 






High quality muon pairs 


1051593 


739947 


Pt > 1 GeV/c, Ip > 1.5ct 


34745 


22359 


mG [2,8] GeV/c 2 ,no J/ij> 


23560 


19254 


Upper limits 


23406 


19129 


P x > 0.04, L k < 0.9 


16750 


13268 


Ip > 4(7 


582 


402 


Az > 0.4 cm 


167 


100 



Table 2. Double muonic b decay selection and surviving 
and /i ± /i ± pairs in the data obtained from cut optimisation. 



In the main measurement method, the Az cut is relaxed 
to 0.2 cm (see Section |73T >. For this cut value, the result- 
ing numbers of opposite-sign and same-sign muons are 225 
and 117, respectively. 

The invariant mass distribution of dimuon events sur- 
viving the bb selection is shown in Figure |2] (the invariant 
mass cut on the J/ip is removed for illustrative purposes). 

7.3 Cross section determination 

In the main measurement method, the number of bb decay 
events in the data (iVft ~) is obtained from a likelihood fit 
to the data of the simulated pr and Ip distributions of sig- 
nal and background events. The selection criteria listed in 
Table [2] are applied, with the exception of the lower limit 
on Az, which is decreased to 0.2 cm. The number of sur- 
viving [i + fi~ pairs becomes 

N+- = 225±15 stat . 

The selection on Az is relaxed because the likelihood 
fit uses more information than that used in the cut opti- 
misation. The likelihood fit is also sensitive to the shapes 
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Figure 2. Dimuon invariant mass after bb selection. Events in 
the J /if> mass region ([2.95, 3.25] GeV/c 2 ), which is highlighted 
in the picture, are mostly due to 6— > J/ip + X decays. Since 
they are used in a different analysis [6|, they are removed. 



of signal and background distributions, while the cut op- 
timisation only uses the number of signal and background 
events surviving the selection. 

The likelihood function is defined as 



L(n s ,n b ) = 



(n a + n b ) n e-^ +n ^ 



nl 



i=l 



n 



S P S + n b P b 
n s + n b 



N+~) selected 



The product index i runs over the n ( 
dimuon events, the exponential term accounts for Poisson 
fluctuations of signal and background, n s (= -W^ - ) and n b 
(= N^T + N+~) are fit parameters representing the num- 
ber of signal and background events in p + p~ data, P s and 
P b are the products of the Ip and pr probability distribu- 
tions of the two muons for signal and background events, 
as obtained from the Monte Carlo simulations ll25ll : 

P a = C 8 • (R s Pbi + Pbx) and 

Pb = C& " (R b P C c + Pco)- 

The signal probability (P s ) is the sum of probabilities 
from decays of class Bl (P m ), B2, B3 and BA (P B x) 
with ratio R s = 1. The background probability (P b ) is ob- 
tained by adding the probability for charm and combina- 
torial events (P cS and P co ) with R b = 0.2. The probabil- 
ity ratios R s and R b are determined through Monte Carlo 
simulations. Constant factors (C s and C b ) are used to nor- 
malise the total probability to unity. 



The result of the likelihood fit performed by minimis- 
ing the quantity —2 log L is shown in Figure [3j The mini- 
mum is obtained when n s = 76 ± 12 and n b = 149 ± 15, 
which corresponds to the bb production Cross section: 

o b i = 16.2 ± 2.5 s tat nb/nucleon. 



S 3 

< 



bb -» + X decays ^ 





















10 



15 



20 



o, - [nb/nucleon] 

bb 



Figure 3. Dependence of — 2AlogL = — 2(logL — logL min ) 
on a h i, where L is the likelihood of bb — ► up, + X decays. An 
increase of one unit on the vertical axis corresponds to a ler vari- 
ation of a h i. 



7.4 Cross-check by a counting method 

In a measurement method based on event counting, N^_~ 

D bb 

is obtained by subtracting all possible background events 
from the number of fi + p~ pairs surviving the bb selection 



listed in Table 0(JV+f): 



bb 



N, 



AT 



Assuming that the reconstruction efficiency of bb de 
cays into p + pT is equal to that into p 
class in each target configuration (e bb j i 

which is true to within a few percent, the bb production 
cross section can be written as 



p , for each decay 

e ±± . = e+f. .), 

bb, 



cr, 



bb 



(N+- - AT±±) Aaj/^BRj/^ - a c5 BR c5 X c 



X, 



i,b 



where the term X bb is given by 
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x bi = y, ^j/^,i A i ■ 

The number of dimuon events surviving the 66 selec- 
tion are N+7 = 167 ± 13 and N±± = 100 ± 10. The sim- 
ulated charm background is N^T = 11.7 ± 0.9, while the 
number of simulated bb decays into same-sign muons are 
A r fc ^ ± = 23 ± 7. The resulting number of 66 decays into 
opposite-sign muons is = 78 ± 23 (the uncertainty 

takes into account the correlation between and A r ,^ ± ). 

bb bb 

The result of the event counting method is the 66 cross 
section a bb = 18.2 ± 5A stat nb/nucleon. 

Compared to the result of the likelihood fit, the sta- 
tistical uncertainty is increased. This is due to the fact 
that the cross section is proportional to the difference be- 
tween two numbers having a large statistical uncertainty 
(Nyr~ oc N*~ — Nj~), while, in the previous method, 
iV^t - is obtained from data only. However, for the 

likelihood fit, the selection criteria are relaxed, which im- 
plies that the selected 66 events contain a larger fraction of 
background. The advantage is that the likelihood fit takes 
into account also the shapes of the signal and background 
distributions. 

The likelihood fit is used to estimate the number of 66 
events in the data, while the event counting method serves 
to estimate the systematic uncertainty associated with the 
measurement method. 

8 Systematic uncertainty 

Many systematic uncertainties on the measurement pre- 
sented in this paper are similar to those reported in refer- 
ence |6], where the reader can find a detailed discussion 
of detector and trigger simulations, models used to simu- 
late 6-flavoured hadron production and decay, production 
of 3/ip mesons, fluctuations of the proton-nucleus interac- 
tion rate, beam position and shape. 

The reconstruction efficiency and the production cross 
section of J ftp mesons are determined assuming no polar- 
isation. The systematic uncertainty associated with J/ip 
production includes the effect of a polarization consistent 
with the limits provided by other experiments I3l"ll32ll33ll . 

The contributions specifically affecting the 66 cross sec- 
tion measurement from double muonic 6 decays are evalu- 



ated. Those which are due to event selection are defined as 
the maximum variation of a bb , divided by \/l2, obtained 
when varying crucial quantities in the likelihood fit, such 
as the impact parameter cut (Ip) and the assumptions on 
signal and background composition (R s and R b ). The Ip 
cut is varied between 2 and 6a, R s in the range [0.5, 2] and 
R b in the range [0, 0.4] (Figure 01). 




2 3 4 5 6 

impact parameter cut [ a] 

Figure 4. Dependence of the bb production cross section on 
the muon impact parameter cut. The band represents the vari- 
ation of cross section due to a change of the expected signal and 
background composition (R s and Rf,) in the ranges [0.5, 2] and 
[0,0.4], Respectively. 

For cut optimisation (Section[L2]), a 66 cross section of 
15 nb/nucleon is assumed. Since a variation of the value 
chosen for a bb might result in a different event selection, 
the systematic uncertainty due to the assumption on a bb is 
included in the effect of cut variation. 

The systematic effect due to the measurement method 
is 6%. This effect is estimated by determining a hb with the 
two methods described in Section|7]on the events surviving 
the 66 selection (Table [2]). 

The systematic effect due to the assignment of a prior- 
ity to the classes of double muonic 6 decays is negligible. 

Assuming that all uncertainties listed in Table [3] are in- 
dependent, the total systematic uncertainty is 18%. 

9 Conclusions 

The 66 production cross section in 920 GeV proton-nucleus 
fixed target collisions has been measured by using double 
muonic 6-flavoured hadron decays. The measurement is 
performed with a likelihood fit of the simulated kinemati- 
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Systematic effect 


Uncertainty 


Detector and trigger simulation 


CO/ 
070 


J/tp production models 


2.5% 


bb production and decay models 


5%) 


b lifetime 


1% 


Proton-nucleus interaction rate 




Beam characteristics 


0.5% 




8.9% 


BR(J ftp —* fi^fi ) 


1.7% 


J ftp nuclear suppression (a) 


3.7% 


ill i~ j." 

J ftp event counting 


n 0(7/ 

0.3% 


Efficiency determination 


o no/ 
2.0/0 


Charge factor (6) 


0.3% 


BR{b -► /i + X) 


3.5% 


BR(c -»■ + X) 


3.2% 


BR(6 -► ccs + X) 


5.5% 


Ip 


5% 


R s and R b ratios 


6% 


Measurement method 


6% 


Total 


18% 



Table 3. List of systematic effects in the measurement of er fe g. 
The first six contributions are evaluated in reference (6) . 



cal distributions of the signal and background events to the 
HERA-B dimuon data. The result is 

a bb = 16.2 ± 2.5 stat ± 2.8 sys nb/nucleon. 

The result is consistent with our previous measurement 
a bb = 14.9 ± 2.2 sta t i 2A sys nb/nucleon |6j, which was 
performed on a statistically independent set of events. Both 
measurements are normalized to the prompt J/ip produc- 
tion cross section. 

The combined result of the two HERA-B measure- 
ments, accounting for the correlation between their system- 
atic uncertainties, is 

a bl = 15.4 ± 1.7 stat ± 1.2™" ± nb/nucleon, 

which is consistent with the latest QCD predictions of Bon- 
ciani et al. [ 1 ] and Kidonakis et al. [ 2 ] based on NLO cal- 
culations and resummation of soft gluons (Figure [5]). 



e 
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" * E789 (b-»J/y + X) 

▼ E771 (bb^nn + X) 

: • HERA-B (b^J/y + X and bb^nn + X) 








1 Kidonakis et al. (2004) 

Bonciani et al. (2002) 







600 700 800 900 1000 

Proton Beam Energy [GeV] 



Figure 5. Cross section for bb production as a function of the 
proton energy in fixed target collisions. The predictions of Bon- 
ciani et al. Q| and Kidonakis et al. |2| are shown. The theo- 
retical uncertainties are obtained by changing the renormalisa- 
tion and factorisation scales and the b mass. The HERA-B re- 
sult, which is based on a combined analysis of b — * 3/tp + X 
and bb — * /i/i + X decays, is consistent with the theoretical pre- 
dictions. The results of the lower energy Fermilab experiments 
(E77 1 and E789 H) are also shown. 
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